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ABSTRACT 




Existence of solutions of two-dimensional "boundary- value 
problems of the type 

= 0 in R, 

u(P) =0 on r x . 


= y(P,u) on r 2 , 

is discussed under conditions subsequently strengthened for 
uniqueness. In general, solutions are shown to lie between a 
maximum possible and a minimum possible, which extremes are 
obtainable as limits of sequences of solutions of certain lin- 
ear problems. Under convexity assumptions on y, the unique 
solution is shown to be obtainable as a maximum over solutions 
of an entire class of linear problems, and is also obtainable 
by Newton’s method. 

Finite difference approximations are shown by a 
Gerschgorin-type analysis to converge to solutions of the above, 
at essentially the rate OCh 1 / 3 ) as h -»• 0, where h is the 
maximum mesh width. The order of convergence is tied to the 
smoothness of y, r^, and r 2 , and more general results are 
actually obtained. 


til t; 

A.P 



The two-dimensional Laplace equation was used for simplic 
ity and for conciseness of treatment. Generalizations are in- 
dicated in the direction of higher dimensions, more general 
elliptic operators, non- homogeneous, non-linear differential 
equations and conditions other than Dirichlet on Tq. 
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CHAPTER I 


INTRODUCTION 


1.1 Object 


Mixed "boundary problems of the type 
V 2 u = 0 in R, 
u = 0 on IN, 


— = r(u) on r 2 , (1.1.3) 

in two dimensions are investigated. R is a bounded two- 
dimensional region with boundary U T 2 . Existence of solu- 
tions is obtained under the mild assumption 

lim sup r(ul < o, (1,1.4) 

I u| -» 00 u — 

which, in effect, permits the replacement of y by a bounded 
function y , having the same values as y within the a-priori 
estimates on u made possible by (1.1.4). Uniqueness is es- 
tablished under the additional, more restrictive hypothesis 

y u <o. (1.1.5) 

The problems of existence, uniqueness, and numerical ap- 
proximations of solutions of Dirichlet problems for elliptic 
equations with solution-dependent source terms, e. g. 
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| V 2 u = f (u) (1.1.6) 

have received a great deal of attention recently and a number of 
results have been obtained, for example, in [1, 4, 9, 12, 13, 


14, 18, 20]. 

Problems like (1.1.1 - 1.1.3) have received some attention, 
for example, in [24] where they are treated from the point of 
view of the variational calculus and (1.1.3) is "weakly" satis- 
fied, but I have not seen any work on the numerical approxima- 
tion to their solutions with estimates of the error in the 
appr oximat i on. 

Techniques similar to those applied to Dirichlet problems 
with a non-linear source term and to linear problems with mixed 
conditions can be applied to yield the desired results for 
problem (l. 1.1 - 1.1.3). Thus the existence proof has a pro- 
cedure somewhat parallel to that of [18] with results of [14] 
and [17] being called upon. Uniqueness is' demonstrated in a 
manner analogous to that in [9] for (1.1.6). The proof that 
the solutions of ( 1.1.1 - 1.1.3) can be approximated arbi- 


trarily closely by solutions of corresponding sets of differ- 
ence equations is a Gerschgorin-type of analysis. The solu- 
tions being permitted discontinuities where Tp intersects 
Tg, the error estimates must follow lines similar to those set 


down by [22] and [21] for Dirichlet problems with non-smooth 


boundaries. The result here is 0 


jpO 2 / ( 2+a) 


convergence where 
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a is the Holder coefficient in the smoothness hypotheses on 
y , It, and f 2 • 

CHAPTER II contains the above-mentioned results on exist- 
ence and uniqueness of solutions u of (1.1.1 - 1.1.3) and, 
under suitable further restrictions on y, a proof of a maxi- 
mum operation for u similar to that established in [13] for 
(1.1.6). CHAPTER III contains the method of finite differ- 
ence approximations and results for them analogous to those 
of CHAPTER IL The error analyses are presented in CHAPTER IV, 
and some possible extensions of this paper in CHAPTER V. The 
results of CHAPTERS II through IV are then summarized in 
CHAPTER VI. Some numerical results are given in CHAPTER VII. 

The remaining portion of this chapter is devoted to two 
problems from the transfer of heat which can be treated by 
proper modifications of the methods given in this paper. A 
special case of the first is used as a numerical example 
illustrating the results. 

1. 2 Motivation 

The first example of a problem which can be put into a 
form suitable for analysis by the methods of this paper arises 
from requirements of rejection of waste heat into space by 
radiation. The heat transfer here is governed by the Stephan- 
Boltzmann law and the mathematical model yields 

V 2 T =0 in R, 


T = % 


on ri. 


( 1 . 2 . 1 ) 
( 1 . 2 . 2 ) 
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_ eo 4 
3n " “ K 


on r 2 . 


(1. 2. 3) 


where T(P) is the temperature at points P of R, e is 
emissivity, K is conductivity, a is the Stefan- Boltzmann con- 
stant and aT/3n is the derivative in the direction of the 
outer normal. 

The second example arises from an attempt to keep the walls 
of a combustion chamber from melting. In many heat transfer 
problems arising in combustion, cryogenics or in those arising 
in combining both technologies, large temperature gradients are 
frequently encountered. The solution of the equations arising 
from the mathematical model assuming constant material prop- 
erties no longer represents the temperature at points in the 
heat- carrying body. 

Assume that part of the boundary is kept at a constant 
temperature while the transfer of heat at the remaining portion 
is governed by Newton's law. Let the conductivity K >0 be 
a function of temperature. Then the equations for the tempera- 
ture are 

div Jk(T) grad t] = 0 in R, (1.2.4) 

T = T c on Hp (1.2.5) 

K(T) = h(Tg - T) on r 2 , (1.2.6) 

where h > 0 is the heat transfer coefficient and, typically, 
Tg, the temperature of the gas, is such that Tg » T c . 
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Under the transformation 


*T 


u = 


K(T) dT 


(1. 2. 7) 


which is 1 - 1, the boundary value problem becomes 

V 2 u = 0 in R, 
u = 0 on Ti, 


f“ = h [ T g “ T < u >] on r 2* 


( 1 . 2 . 8 ) 
(1. 2.9) 

( 1 . 2 . 10 ) 


Thus (1.2.4 - 1.2.6) can also be placed (by (1.2.7)) into a form 
suitable for treatment by the methods of this paper. 

The advantage of having a problem of this type in this 
latter form becomes apparent upon examining the system of finite 
difference equations. The system corresponding to (1.2.4 - 
1. 2. 6) has most of its equations non-linear while that corre- 
sponding to (1.2.8 - 1.2.10) has most of them linear. 


CHAPTER II 


SOLUTIONS OF THE BOUND AEQ" VALUE PROBLEMS 
FOR THE PARTIAL DIFFERENTIAL EQUATIONS 


2. 1 Existence 



Consider 

V 2 u =0 in R, 

(2.1.1) 


u(P) =0 on F 1 , 

(2.1. 2) 


~ = r(P,u) on r 2 , 

(2.1.3) 


where R is a hounded, two-dimensional region with boundary 
I'd U Fg, consisting of a finite number of smooth arcs. Here 
Fj denotes the closure of F]_. As was done for Dirichlet 
problems with non-linear source terms by Levinson [14] and 
Parter [18], r(P>u) is to be replaceable by a bounded function 


Tq by requiring that 


lim sup r(P>u) < 0 
| u| -*• oo u — 


(2.1. 4) 


Before proceeding to smoothness requirements on y and on 


the boundary of R, it is necessary to make some definitions. 

A point P of the boundary U is said to have a 
barrier function there, if there exists a superharmonic function 
cDp, continuous and single valued in R, such that 

(a) cop(P) = 0, 

(b) o>p(Q) * 0, for QeR - {P}. 
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T 1 is said to be smooth if it can be covered by a finite 
number of circles in each of which, one of the co-ordinates can 
be expressed as a function of the other, the "arc parameter". 
These functions are required to have Holder continuous second 
derivatives . 

Some smoothness in the form of Holder conditions on deriva- 
tives of T with respect to arc parameters and with respect to 
the dependent variable is also required. Namely, let there be 
a number a such that 0 < a < 1, and five positive functions 
K(M) ,Kq(M) ,K p(M) ,Kg(M) , and K^(M) such that for all P,Pp, and 
P^er^ and all u, Up, and u^ whose absolute value is bounded 
above by M, 


I r p( p l' u ) " 

r P (P 2 ,u)| 

£ 

K(M)|P 1 - 

P ^ 

21 > 

(2, 

.1.5) 

|r(P,up) • 

- r(P,u 2 )| 

£ 

K 0 (M>h 1; 

- u 2l> 

( 2 , 

.1.6) 

l r u (P l’ u) - 

r u ( p 2 ,u)| 

£ 

K 1 (M)|P 1 

“ *2^ 

(2, 

.1.7) 

l r u (P,u l) - 

r u (P ’ u 2> 

£ 

*2< M >K 

' u 2 1 ’ 

(2, 

.1.8) 

|Tp(P,Up) - 

Vp( p ? u 2 ) | 

£ 

K 3 (M)| u i 

‘ u 2h 

(2, 

.1.9) 


where Tp denotes the partial derivative with respect to arc 
parameter . 

Some additional restrictions of dR are necessary to in- 
sure existence, uniqueness, and some degree of smoothness of 
solutions of some auxiliary problems of this section. Namely, 

(i) A barrier function must exist at each point of Pq U Tg . 

(ii) The tangents to Tp and Tg at the points of Tp H 

must not meet at either of the two angles 0 or n . 
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(iii) must "be smooth. 

(iv) is smooth. 

Although used in this section and in section 4.1, this last 
hypothesis can he dispensed with, under modifications to he in- 
dicated. 

It will now he shown that the search for a solution of 
(2.1.1 - 2.1.3) may he carried out in a hounded portion of the 
(P,u) - space. 

Condition (2.1.4) implies the existence of a positive func- 
tion t) ( u) which is decreasing for large u, such that. 


lim t|(u) = 0, 

U-*oo 


( 2 . 1 . 10 ) 


which satisfies, for large u, 

iM. < T](u). (2.1.11) 

u “ 

Furthermore, Levinson showed that t) could he modified, with- 
out destroying these properties, so that ut)(u) + ti is in- 
creasing for large u ( see [14] lemma 3. 1) . 

To derive the a priori estimates, cp is postulated as a 
solution of (2.1.3). The proof, that cp has a hound depending 
only on R and y, depends on the existence of a positive func- 
tion f; for which upper and lower estimates are available, and 
the definition of a new function = cp/1 which is shown to he 
hounded, in two steps: The first is the proof that any positive 

maximum or negative minimum must occur on the boundary Tg. The 
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second step is to show that these extremes have a hound depend- 
ing on 


By Miranda [17], there is a function (; 

continuous in R 

which satisfies 


V 2 £ = -1 in R, 

(2. 1. 12) 

•s 

IC H 

CJ 

o 

1 — 1 
II 

(2.1.13) 

= 1 on r 2 . 

(2.1.14) 


Greenspan [11] and Batschelet [3] use this function to prove 
convergence of finite difference approximations to solutions of 
linear mixed problems under enough conditions that {; will have 
four continuous derivatives. Here, we are only guaranteed that 
the first and second derivatives satisfy Holder conditions in 
compact subsets of R - O T^. 

Let cp be a solution of (2.1.1 - 2.1.3) and define by 

cp = (Mr. (2.1.15) 

Assume that has a positive maximum at an interior 
point P of R. Since at P 

V 2 cp = 0, 

i. e. 

♦ v2? + 2 (H^ + |^) + sv2 * = 0 ' 

we have 

v 2 ^ = 

= i 
i 

> 0 


( 2. 1. 16) 
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at this point. Thus ^ must take on any positive maximum on 
the boundary. Since \|f is zero on Tp this positive maximum 
is possible only on Tg. A contradiction is also reached for 
a postulated interior negative minimum. 

Let M-j_ > 0 be such that for u > M^, T](u) is decreasing 

and 

r 2 

and 


r(p>u) 

u 


< T)(u) • 


(2.1.19) 


Now suppose Mg, the positive maximum of t at PeTg is such 
that 


^ >M r 

Then at P, by (2.1.3), (2.1.14), and (2.1.15), 



< 


( 2 . 1 . 20 ) 


( 2 . 1 . 21 ) 


( 2 . 1 . 22 ) 


which contradicts the maximality of \Jr at P. Thus (2.1.20) is 
false, i.e. 


Mg <M r 


(2.1.23) 



11 


Similarly a negative minimum for if on Tg has smaller magni- 
tude than M-j_. 

So, by (2.1.15), any solution cp of (2.1. 1 - 2.1.3) has 
the a-priori bound 

|9| <% = m | X ? • Mi (2.1.24) 

which depends only on R and rj. This establishes the fol- 
lowing theorem: 

Theorem 2. 1. 1 . - If cp is a solution of (2.1.1 - 2.1.3), 

then cp has the a priori estimate 

< m * x C • M-, 

R 1 



where £ is the solution of (2.1.12 - 2.1.14) and is so 

large that r)(u) is decreasing for u > and (2.1.18) is 
satisfied. 

This paves the way for the replacement of y by a bounded 
function Yq. In fact, let M 4 > Mg + 1 and 


r(P> u) 

t(p,m 4 ) if r(P,M 4 ) <0 n 

Yq(P,u) = /max[ 0, -)r(P,M 4 ) - u + M 4 ] otherwiseju > M 4 (2.1.26) 
r(P, -M 4 ) if r(P> -M 4 ) >0 

max[ 0, y( P, -M 4 ) - u - M 4 ] otherwiseju < -M 4 (2.1.27) 

The hypotheses satisfied by y are also satisfied by Tq. The 
easiest to verify are (2.1.5), (2.1.6). A little more work is 


u i < 


(2. 1. 25) 


required to verify ( 2 . 1 . 11 ), and this is now carried out. 
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For large u but where |u| < M 4 , (2.1.11) is immediate 
for Tq from (2.1.25). Let u > M 4 . Then since ut](u) + u is 
increasing > 

M 4 T](M 4 ) + M 4 < uri(u) + u, (2.1.28) 

so that 

r 0 (P,u) t(p,m 4 ) m 4 

u = u 1 + u 

1+ ^ 

- u U 

hy (2.1.19). Therefore hy (2.1.28) it follows that 

,u) 

<ri(u) (2.1.29) 

if u > M 4 . 

The argument is similar for u < -M 4 . Thus Tq satisfies 
(2.1.11). By examining the several cases which arise it can be 
seen that the transitions at 

u = t(p,-m 4 ) - m 4 , (if t(p,-m 4 ) < 0 ), 
u = -m 4 , 
u = M 4 , 

u = t(p,m 4 ) + m 4 (if t(p,m 4 ) > 0 ), 

can be made smooth enough so that, besides (2.1.5) and (2.1.6), 
Tq satisfies the inequalities (2.1.7 - 2.1.9) without violat- 
ing (2. 1.11). Since Tq does not differ from r for u < ^3 < 
M 4 , cp is a solution of (2.1.1 - 2.1.3) if and only if it is a 
solution of this set of equations with y replaced by Tq. 
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The point of this latter result is that we may proceed to 
investigating existence of the solution of (2.1.1 - 2.1.3) with 
X having one .more (crucial) property: There is an N > 0 such 

that for all PeP 2 and | u| < «>, 

|r(P,u)| < N. (2.1.30) 

Let the distance between two sets be defined as usual: 


d(A,B) = inf [P - Q| , 

P e A 
Q e B 

In what follows, the various functions z, say, defined by 
solving certain boundary value problems, exist by the results 
of Miranda [17]. Here, if the conditions on r 2 are written as 

= $(?)> (2.1.32) 

then the first derivative of $ with respect to the axe param- 
eter needs to satisfy a Holder condition on compact subsets. 
r 2 of r 2 , where d(? 2 , r^r 2 ) >0, and $ to be bounded on 
such sets. The results for z in R then are that z is con- 
tinuous in R and its first and second partials satisfy Holder 
conditions in compact subsets of R - O lb,. 

Let Uq be defined for N* > N by 


V 2 Uq = 0 


u 0 = ° 


in R, 
on r 1# 


duo 

dn 


= N‘ 


Then Uq >0 in R. Let 


(2.1.33) 

(2.1.34) 


on r 2 . 


(2.1.35) 
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and 


Then 


for 


U QM = ma £ * 0 (P) 
PeP 2 


k = k o^ u cm) 


r(u) - r(w) - k(u - w) >o. 


” U CM < u < v < U CM* 

Let u m+1 (P) "be defined in R by 

V2u m+1 = 0 in R > 
u m+l = 0 on 

-^T + k Vl = *(%) + ku m on r 2 ,m>0 
Examine the case m = 0 on 

du- 


+ ku x = r(u Q ) + ku 0 


< N’ + kuQ 
= SU ° + ku 

“SS 0- 


Therefore, 


Also, on r 2 


U 1 ^ u 0* 


_ = -k(u x - u Q ) + r(n 0 ) 

> r( uq) 

> -r 

= _ Su o 

Sn 


(2.1.36) 

(2.1.37) 

(2.1.38) 

(2.1.39) 

(2.1. 40) 
(2.1.41) 

(2.1. 42) 

(2.1.43) 
(2.1. 44) 


(2.1. 45) 
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which implies 


Up > - u 0 . 

Having proved for m = 1 that 


- u 0 < % ^ %-l ^ u 0 


(2.1. 46) 


(2.1.47) 


proceed inductively to the case m + 1: On Tg 

Ih ( %+l - V + k <Vl - V= r(u ^ - ^K-l) + k (% - %-i) 

< 0, (2.1.48) 

so that 

Um+l < % (2.1.49) 

and 

“ir = _k(u to + i ■ V + r( V 

> r(%) 

> -r 

= (2.1.50) 

Sn 

implying (2-1-47) for all m. 

Thus, (u^) is a uniformly hounded, monotonic sequence of 
functions- They are, in addition, harmonic, so that by Harnack’s 
second theorem (Petrovsky [19] p. 178) convergence is uniform in 
every compact subset of R and the limit is harmonic. Since 

]?i U r 2 is assumed to have a barrier function at each point, 
the limit 


= lim ^ 
m-*» 


u 


(2.1.51) 
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is continuous in R U and takes the value 0 on r^. 
(Petrovsky [19] p. 184). Using the results of Miranda, it is 
possible to say more. 

Let S > 0 and R§ be the compact subset of R - 
such that d(R & , Pi F 2 ) = 5: 

R 6 = {PeR|d(P,r i n r 2 ) > 6). (2.1.52) 

Recalling (2.1.42), let 

Vi( p ) - r[P,V p )l + ku m( p ) (2:1.53) 

for m > 0. Then by (2.1.5), (2.1.7 - 2.1.9), the sequence of 
first derivatives of the $ with respect to arc parameter 
satisfies a Holder condition uniformly in Rg O Tg. By 
Theorem 6,1 of Miranda [17] the second partial derivatives 
satisfy uniform Holder conditions in Rg and, therefore, the 
sequence of second partial derivatives of form an equi- 

continuous set. There is, then, a subsequence of second deriva- 
tives which converges uniformly, so that the boundary conditions 
on Pg are satisfied by the limit u. Furthermore, the first 
and second derivatives of u are Holder continuous on compact 
subsets of R - H Tg. 

These last results can be summarized as: 

Theorem 2. 1. 2 . - If R is a bounded region of the xy-plane 
with boundary U r 2 > satisfying the smoothness conditions (i) 
through (iv) on page 7, and y(P, u) is a function defined for 
all PeFg and all finite u, satisfying the smoothness condi- 
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tions (2.1.5) through (2.1.9) on page 7 and the fundamental 
condition (2.1.4) on page 6, then there exist solutions of 
(2.1.1 - 2.1.3) continuous in R and having Holder continuous 
first and second derivatives in compact subsets of R - 

The bases of the existence theorem 6,1 of Miranda are 
some Schauder-like estimates on spaces of functions which, 
with some of its derivatives, are permitted various orders of 
discontinuities as points approach Pi Tg. If it is desired 
to relax the smoothness of (to dispense with (iv) on page 8), 
Miranda has other estimates which become more nearly like the 
usual Schauder estimates (see Courant and Hilbert [9] p. 331 ff.). 
The above arguments are modified, chiefly in the definition of 
Rg, namely, 

R 5 = {PeR|d(P,r 1 ) >6). (2.1.54) 

This gives rise to a modification of Theorem 2.1.2 where the 
compact subsets in the conclusion are those of R - r^. 

The proof of this theorem could also call upon results of 
Agmon, Douglis, and Nirenberg [2], especially on Theorem 12.2, 
but more care would have been necessary in defining and 

cp^, since the boundary operators need to be smooth on the 
entire boundary. 

The solution u will be seen in what follows to be a 
maximum possible solution to (2.1.1 - 2.1.3). In a similar 
fashion a minimum possible solution v can also be obtained. 
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Briefly, (omitting the details as they are parallel to those 
for u) define v Q by 

V 2 v Q = 0 in R, (2.1.55) 

v 0 = 0 on (2.1.56) 

dv 0 

= - N * on r ?- (2.1.57) 

dn c — 

Define k analogously to (2.1.37) and (v m ) by (2.1.46 — 2.1.42). 

Theh 


v = lim v m . 

rrnco 

Let Z be any solution of (2.1.1 -2.1.3). On T 

I s = r(z) 

Sn 

< N‘ 

Letting w = u Q - Z, by (2.1.60), since 6u^n = K», 

> 0 

bn 


(2.1.58) 

(2.1.59) 

(2.1.60) 

(2.1.61) 


on r^ t so that 


Z < Uq in R. 


(2.1.62) 


To show that Z is bounded above by all iterates u m another 
inequality implied by (2.1.6), (2.1.37), and (2.1.39) is needed, 
namely, if -u^ < u < w < u QM then 

r(w) - r(u) - k(h - w) >0. (2.1.63) 

Let w = Up - Z. Ey ( 2. 1. 43) , on r 2 

6up 

— = r(u n ) + k(u 0 - up) 
on w 

= r(u 0 ) + k(uQ - Z) + k(Z - Up) 

= r(u 0 ) + k(uQ - Z) - kw. 


(2.1.64) 
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Subtract (2.1.59) from (2.1.64) and obtain 

~ + kw = r(u 0 ) - r(z) - k(z - U 0 ) 

>0 (2.1.65) 

by (2.1.62) and (2.1.63). Therefore 

Z < u-l. (2.1.66) 

Replacing u Q by u m and u-j_ by u fil+1 completes the induc- 
tion. That is, m > 0 implies 

2 < Ujj (2.1.67) 

for any .solution Z of (2.1.1 - 2.1.3). The result is 

Z < u. (2.1.68) 

Carrying out a similar argument with v obtain finally that 

v < Z < u, (2.1.69) 

for any Z solving (2.1.1 - 2.1.3). 

These final results may be summarized as 

Theorem 2.1.5 . - Under the hypotheses of Theorem 2. 1. 2, if 
the solution u is defined as in its proof, namely, by 
(2.1.51), and the solution v by (2.1.58), then v<u and 
if z is any solution of (2.1.1 - 2.1.3), then 

V < z < u. 

This result is like that obtained in Courant and 
Hilbert [9], pp. 369-372, for Dirichlet problems for 

V 2 u = f(P,u) in R, 

and later by Parter [18] under more general conditions on f 
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and the boundary of R, still however , for Dirichlet problems 
for the non-linear Poisson equation directly above. 

2.2 Uniqueness. 

In addition to conditions (2.1.4 - 2.1.9), imposed to 

guarantee the existence of solutions to problem (2.1.1 - 

2.1.3), further restrictions must be placed on y to guarantee 

the uniqueness of a solution . For an example of a boundary 

value problem satisfying all conditions of section 2.1 but 

possessing infinitely many solutions, consider the two-point 

boundary- value problem , 

.2 

SJi =0, 0 < x < 1, (2.2.1) 

dx 2 

u = 0, x = 0, (2.2.2) 

§£ = r(u), x = 1, (2.2.3) 


where 







u < -1 







r(u) =< 


-1 < u < 

1 


(2.2.4) 





w 1 ' 

u > 1. 




Here one may 

take 

K(M) s ; 

1, any a, 0 < a 

< 1, 

K 0 (M) 

= 1 and 

T) ( U ) = 

l/u. 

The 

solution 

i 

of (2.1.12), 

(2.1 

• 13), 

(2.1.14), 

where 

Fi is 

> (0) 

is 









£(x) 

= 1 

+ 2x - | x 2 . 



(2.2.5) 


so that 

max £=2.5 y (2.2.6) 

R 
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and so 


M x = 2.5, 

and the a priori hound (2.1.24) for solutions of (2.2.1), (2.2.2), 
(2.2.3) reads 

|cp| < M 3 = 6. 25. ( 2 . 2. 7) 


But there are many solutions of this problem, seme of which are 
given by 

u = ax (2. 2. 8) 

for any | a| < 1. The constant N of (2.1.30) may be taken to 
be N = 1. The function Uq of (2.1.33), (2.1.34), (2.1.35) is 

Uq(x) = x, (2. 2.9) 

so that Uq^ = 1, k = 1, and therefore for all m > 0, 

^(x) = x. (2. 2. 10) 

Therefore, the limit is 

u(x) = x, (2. 2. 11) 


and similarly 

v(x) = -x. (2.2.12) 

So if z is any solution of (2.2.1), (2.2.2.), (2.2.3), 

-x < z(x) < x, 0 < x < 1. (2. 2. 13) 

In order to obtain uniqueness it is further required that 
T satisfies 

r u (P,u) < 0. (2.2.14) 

Other similar conditions also can be used to deduce uniqueness. 
For example, see the papers of Ablow and Perry [1], Kalaba [13], 
and Pohozaev [ 20] . 



22 


First, consider the linear problem 

V 2 w = 0 in R, (2.2.15) 

w = 0 on Fp, (2.2.16) 

§5 = -c(P)w on r 2 , (2.2.17) 

where c > 0. That any postulated positive maximum or negative 
minimum for w must he assumed on r 2 follows from (2.2.15) 
and (2.2.16). If c > 0, then either of these postulates con- 
tradicts (2.2.17) at the point P* s Pg where the postulated 
maximum (minimum) occurs. If at P*, c = 0, then consider \|r 
defined by 

w = £* (2. 2.18) 

where £ is defined by (2.1.12- 2.1.14) so that as before, 

(sec. 2.1) if i|r has a positive maximum (negative minimum) it 
occurs on Tg. But by (2.2.18), (2.1.14), (2.2.17), 

1^- = -%■ at P* (2. 2.19) 

3n (; 

= P/l*. (2.2.20) 

For either postulate, cty /Sn will have the wrong sign. 

Therefore, the following lemma is obtained: 

Lemma 2. 2. 1 . - The linear problem (2.2.15 - 2.2.17) with 
c(P) >0, has only the one solution w = 0. 

Next, suppose the non-linear problem (2.1.1 - 2.1.3) has 
two solutions u and v and define 

w = u - v. 


( 2 . 2 . 21 ) 
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which satisfies 

V 2 w = 0 in R, (2.2.22) 

w(P) =0 on F 1 , (2. 2. 23) 

= r(u) - r(v) on r 2 . (2.2.24) 

By the mean- value theorem, (2.2.24) becomes 

= -c(P)w on r 2 , (2.2.25) 

where 

c(P) = -r u (P,cp) >0, (2.2.26) 

and cp is some number between u(P) and v(P). By Lemma 2.2.1, 
problem (2.2.22), (2.2.23), (2.2.25) has only the solution 
w = 0, i. e. 

u = v. (2. 2. 27) 

It is possible to get uniqueness even if y is not dif- 
ferentiable, by replacing (2.2.14) by 

if u < v then r(v) < r(u) . (2.2.28) 

With this and (2.1.6), if u(p) ^ v(P), the boundary condition 
(2.2.25) still holds with 

c(P) = GK 0 (M) >0 (2. 2. 29) 

for some y such that 0 < 6 < 1. 

These results can be stated as 

Theorem 2. 2. 1 . - Under the conditions of Theorem 2. 1. 2, if 
r(P,u) is non-increasing as a function of u, in particular if 
< 0, then there is only one solution of (2.1.1 - 2.1.3). 
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The proof is similar in spirit to that of Courant and 
Hilbert [9], pp. 320 - 324. 

2.3 Newton 1 s Method and the Maximum Operation . 

By strengthening the hypotheses on y, results may be ob- 
tained which parallel those of Kalaba [ 13] where solutions of 
Dirichlet problems for quasi-linear differential equations, for 
example 

V 2 u = f(P,u) in R, (2.3.1) 

were obtained by Newton's method as a maximum of solutions of a 
class of linear equations. 

For problem (2.1.1 - 2.1.3), the conditions of Theorem 
(2.1.2) are to be supplemented by y u < 0 and 

T(u) = max[y(v) + (u - v)r'(v)], (2.3.2) 

v 

where the prime indicates the partial derivative with respect to 
u, i.e., y is convex in u. That y satisfies the usual in- 
equality for convexity can be seen by adding the following two 
inequalities resulting from (2.3.2): Let w = (l/2)(u + v). 

Then 

r(u) > yjjj- (u + v)J + (u - w)r'(w), 

and 

r(v) > rjj- (u + v)] + (v - w)r'(w). 

A sufficient conditions for convexity is that y be twice 
continuously differentiable in u and y uu >0. 
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The property of u being a maximum is altered to that of 
being a minimum if y is concave, i.e. , if 


r(u) = min y(v) + (u - v)r'(v) . (2.3.3) 

v L 


For the concave case, however, it is not necessary to 
separately require (2.1.4) since it is then a consequence of 
(2.1.14) and (2.3.3). 

Newton’s method applied to the problem of (2.1.1 - 2.1.3) 
is carried out as follows. Let uq be any initial approxima- 
tion, good or poor, to u, the unique solution of (2.1.1 - 
2.1.3) guaranteed by sections 2.1 and 2.2. Consider the 
sequence of approximations {i^} defined by the linear problems, 
m > 1 

V 2 u m =0 in R, (2. 3. 4) 

Ujjj = 0 on r 1} (2. 3. 5) 

chu 

-v- = r( Um.q) + r'(u m _ 1 )( u m - u^) on IV,. (2.3.6) 
On 

To compare u with let w = u - ^ and note on r 2 that 

^ = r(u) - [W^.i) + r'(%_i)(u m - u m _ 1 ^ 

^^u” 1 - 1 ) + r'(u m „i)(u - U m _ 1 ) - ^(Ujn.i) + r'K.pXu^u^j 
by (2.3.2). That is 

V^w =0 in R, 

on T lf 


w = 0 


(2.3.7) 

(2.3.8) 
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‘ T ’( u m-l) w > 0 on r 2* (2.3.9) 

Therefore, because of (2.2.14), w > 0, i.e., for all m > 1 
Ujjj < u in R (2.3.10) 

independently of the choice of Uq. 

To compare with m > 1 note that on r 2 

Sil. 

3JT = ^"tn-l^ + ’ Mm-lJT'Cum-i) 

<r(V (2.3.11) 

by (2.3.2) and 

du 

= + (u m+l - % )TT’(%)* (2.3.12) 

Subtracting (2.3.11) from (2.3.12) obtain for w = - u 

that 

- r’tujv >0 on r 2 

whereas in R, (2.3.7) holds and on F^, (2.3.8) holds, so that 
again w > 0. That is 

u m+l - % in R * (2.3.13) 

With this and (2.3.10), is a uniformly bounded, monotonic 
sequence of harmonic functions and, as in section 2.1, its limit 
solves (2.1.1 - 2.1.3). By Theorem 2.2.1 this is the same func- 
tion defined by the limit of the sequence defined in section 2.1. 
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Since this problem is equivalent to (2.3.4 - 2.3.6) with 
Uq = u (so that u^_ = u) , (2.3.10) implies that 

u = max juj (2.3.4 - 2.3.6) with m = lj . (2.3.14) 

If y were instead concave, a similar result with "max’ 1 re- 
placed by "min" in (2.3.14) would hold for u. Here, of 
course, (u m ) is a non- increasing sequence. 

Summarizing, 

Theorem 2. 5. 1 . - Under the conditions of Theorem 2. 2. 1, 
namely smoothness and 

r u - °’ 

if T is convex, then the solution of (2.1.1 - 2.1.3) may he 
obtained by Newton’s method and the maximum operation (2.3.14) 
is valid. 


CHAPTER III 


APPROXIMATE SOLUTIONS 

5.1 Derivation of the Difference Equations . 

In this section, difference approximations are made which 
lead to systems of algebraic equations whose solution is shown 
in Chapter IV, to differ from the solution of the original 
boundary value problem by an amount tending to zero with h at 
the rate h!/3. 

Higher order approximations are discussed in 

Chapter V. 

In what follows it will be convenient to consult figure 1 
where R^ is displayed. Here has its points circled, 

Pj has its points with asterisks superimposed, and R^ has 
its points merely emphasized. Dashed lines on this figure in- 
dicate sections which are singled out for further treatment in 
figures 2 and 3. 

To derive the difference equations, a square net of mesh 
width h is superimposed on the bounded region R and a sys- 
tem of equations is obtained, one equation at each node of the 
net. These equations differ, depending upon the relative posi- 
tion in the net and three types must be distinguished according 
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as the node lies in R, near F-j_ or on r 2 - 

All intersections of U with grid lines will he 
considered along with intersections, interior to R, of hori- 
zontal grid lines with vertical grid lines, as nodes of the 
net. 

The totality of all nodes is denoted hy R^. Each point 
of R^ will have its corresponding difference equations (de- 
generating on r i to a mere statement of the functional value 
at that point). The notation used here is like that of various 
papers of Bramble and Hubhard (see, e. g. , [7]), Let (r=T^ur2) 

r h = rnR h (3.1.1) 

r* = Jp PeR h - r h but P has l 

I a neighbor QeT h J (3.1. 2) 

R h = \ " (r h U (3.1.3) 

The sets and r* are further decomposed into 

r hl = r h n Tq, (3.1.4) 

and 

P h2 = P h O r 2‘ f 3 ' 1 ' 5 ) 

* * 

r hl and are defined analogously except here the repre- 

sentation 

J 'h = I hi'J r h2 (3.1.6) 

would not necessarily be a disjoint union unless it is re- 
quired that those points P which have neighbors in both 
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Figure 2. - Points used in interior 
approximations . 



Figure 3. - Points used in ■boundary- 
approximations . 
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r hi and 


r h2 “be placed in P^, thus making 


r hi n - 0- 


(3.1.7) 


Figure 2, with hq < h, exhibits 5 points singled out by dashed 
lines in figure 1, and labelled 0, 1, 2 , 3, and 4 in this new 
figure. Here, for example, 0 corresponds to the point Pq = 
(x 0 ,y 0 ), 1 to the point P 1 = (x Q + h-^yo) etc. Let 
Pq e Rh U r.. Then either through Taylor's series (see 
Forsythe and Wasow [ 10] pp. 179 ff . ) or through an integration 
method (Varga [26] pp. 181 ff . ) one obtains at the point Pq 
the equation 


2ui 


2u- 


V h l h 3 


h 2 h 4 J 


u~ - 


0 " hjThi + h 3 ) “ h 2 (h 2 + h 4 ) 


2u7 


2u/ 


h 3 ( h 4 + h 3 ) h 4 (h 2 + h 4 ) 


= 0, 


(3.1.8) 


corresponding to (2.1.1). For P 0 € T l±> one or more values 
(say u 4 ) are known to be zero and (3.1.8) reads 


2u n 


2u^ 


2u, 


= 0. 


\h 1 h 3 + h 2 h 4 J U 0 " hpCh-L + h 3 ) - h 2 (h 2 + h 4 ) ' h 3 (h 3 + hq) 

(3.1.9) 

For Pq e r^q the equation is merely 


u 0 = 0. 


For Pq e the equation is the usual 

(4u 0 - uq - u 2 - u 3 - u 4 ) 


= 0 


(3.1. 10) 


(3.1.11) 


h 


h2 


leaving only points in P- 


to consider. 
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Following Greenspan [ 11] , either the "boundary normal ex- 
tended inwardly from Pq meets a grid line obliquely at 
another node labelled P^, or between two nodes. In the first 
case, the outwardly directed normal derivative is replaced 
simply by 

u o - U 1 

d > 

where d = 0(h) as h -*■ 0 is the distance between Pq and P]_. 

In the second case, consider figure 3. The normal derivative 

at Pq is replaced by 

u 0 d i ^ d 2 

d " d(d x + d 2 ) U 2 " d(d x + d g ) U 1‘ 

The equations at these points are thus in one of two forms, 

Uf> - Ut 

d - = r(p 0 ,u 0 ) (3.1.12) 


or 


u, 


0 


d l U 2 


d 2 U l 


d(d x + d 2 ) " d(d^ + d 2 ) r ( p 0' u C))' 


(3.1.13) 


5. 2. Solution of the Difference Equations . 

The system of equations obtained by combining equations of 
type (3.1.8 - 3.1.13) written at all points of R^ where they 
apply may be considered in the form 

Au = d(u) (3.2.1) 

where the matrix A is such that det A ^ 0 and in fact, A is 


an M-matrix so that 
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A -1 > 0. (3. 2. 2) 

This follows from a lemma due to Varga (see Roudebush [22] 
lemma 2, p. 11) since < 0 for i ^ j, aj_^ > 0, and A is 

reducible but has normal reduced form (Varga [26], p. 46) with 
irreducibly diagonally dominant submatrices along its diagonal. 

Let N be defined by ( 2. 1. 30) and consider Uq defined 


■by 


AUq = b (3. 2. 3) 

where bj_ is the same as the i^ component of d(u) in (3.2.1) 
unless Pi 6 ^^’ At these points, (3.1.12) (analogously 
(3.1.13)) is replaced by 


u 00 " U 01 
d 


N + 1. 


(3. 2. 4) 


Then Uq _> 0. Let 
where ||uq|| denotes 


k = K 0 (||u 0 ||), 


max| uq, i | • 

i 


(3.2.5) 


Then 

r(u) - r(v) - k(u - w) > 0 (3.2.6) 

for 

- \\}±o\\ < u < w < Ihioll- (3.2.6a) 

Define u , n , m > 0 by 
— m+1 — 


A iVl = 


(3. 2. 7) 
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where the i th component of b^) is the same as that of 
d(u) of (3.2.1) unless P^r^. Here 

i = £(2m) i + k K,i - %fl,i) 

= rCu^) + k( %#1 - (3.2.8) 

This system of equations can also he written as 

\+l = IK) + (3.2.9) 

where 

K i = °> p i i r h2 

*i = ku mi' P i € r h2> (3.2.10) 

and A differs from A only in the corresponding rows by 
having k added to the (already positive) diagonal element in 
those rows, so that 

0 < A” 1 < A” 1 . (3.2.11) 

Let Pq e r h2 . Then for m = 0, (3.2.9) reads 

0 ” U 1 1 

“^“d — — + ku i, o = r(u 0 ^ 0 ) + ku 0 ^ 0 

< N + ku Q? o+l 

r U 0.0 ~ U 0.1 + ku 

d Ku 0,0. (3.2.12) 

A like inequality holds if the equation at P Q were of the 
form (3.1.13). At other points of R^, 

X a ij U m+l,j = = °* (3.2.13) 

If (3.2.12) and (3.2.13) are combined, one obtains 
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Aui < Au q , (3.2.14) 

from which follows 

Uf < Uq. (3.2.15) 

If (3.2.9) is again examined in the light of (3.2.15), one ob- 
tains 

U 1 0 - U 1 1 

d ~ = ^ u O,o) + k ( u 0,0 - u l,o) 

> r( u 0>0 ) 

> -N - 1 

• ~ ^° ' U °'t (3.2.16) 

hy ( 2. 1. 30) . Therefore, 

Au-^ > -Auq, (3. 2.17) 

so that 

"Ho ^ Hi* (3. 2 * 18 ) 

If (3.2.15) and (3.2.18) is applied to (3.2.6) one obtains 

A(ug - u x ) < 0 (3. 2.19) 


or 


2*2 ^ Hi- 


(3. 2. 20) 

This and (2.1.30) implies (as (3.2.17) followed) that 

Au 2 > -Au 0 (3.2.21) 


or 


-u 0 < n 2 - (3.2.22) 

Repeating the steps leading from (3.2.15) and (3.2.18) to 
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(3.2.20) and (3.2.22) with the subscripts 1 and 2 replaced by m 
and m + i, m > 1, completes the induction, i. e. 

-*0 ^ 2Vl < *m ^ * • * < *2 £ Hi £ 2Jo* (3> 2 * 23) 

Therefore, the components of (u^) form a bounded monotonic 
sequence of real numbers and so have limits. Since only a 
finite number of arithmetic operations are involved in each equa- 
tion of (3.2.7), the limit vector 

u = lim (3.2. 24) 

ITHoo 

satisfies (3.2.1) by virtue of (3.2.8). 

Similarly defining v^ by 

Avq = b (3. 2. 25) 


where instead of (3.2.4) 

v °° _ - vq, J- _ _ ! (3.2.26) 

is satisfied, and v m hy (3.2.7) for m > 1, obtain a monotoni- 
cally increasing sequence whose limit 


v - lim 

ITHoo 


(3.2.27) 


also satisfies (3.2.1). 

It is now shown that any solution of (3.2.1) is bounded from 
above by u and from below by v. To this end, let z be any 
solution of (3.2.1). Again comparisons of difference expressions 
are carried out only on T-^y the only place they might differ. 
Subtracting the latter of 
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= r(* 0 ) 

< W + 1 

from (3.2.4) obtain for w = Uq - z 


W 0 “ W 1 


> 0 . 


(3. 2.28) 
(3. 2. 29) 


(3.2.30) 


Therefore 

z < Uq. (3. 2. 31) 

One more inequality beside (3.2.6) is implied by (3.2.5), 
(3.2.6a) and (2.1.6), namely, if - ||uq|| < u < w < ||uq||, then 

y(w) - y(u) - k(u - w) > 0. (3.2.32) 

Wow let w = Ui - z. By (3.2.8) 

u 10 " U 11 

d = r( u 00 ) + k(u 00 - u 10 ) 

= r(u oo) + k < u oo ■ z o) + k ( z 0 - U 10^ 

= r(u 00 ) + k(u 00 - Zq) - kw Q . (3.2.33) 

Subtracting (3.2.28) from this one obtains 

Wq - W-j^ 

— g — + kw 0 = r(uoo) - r( z 0 ) - k(z 0 - Uqq) 

> 0 (3. 2.34) 


by (3.2.31) and (3.2.32). But this implies 

Aw > 0, 

i. e. 


(3.2.35) 


z_ — Vh • (3. 2. 36) 

Induction is completed in the same way as before for (3.2.23) 
and the conclusion is that 
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z < u. (3. 2. 37) 

Similarly, 

v < z. (3. 2. 38) 

That is, any solution of (3.2.1) is in the "interval" [v,u]. 
Therefore, the following theorem is established: 

Theorem 5. 2. 1 . - Under the conditions of Theorem 2. 1. 3, 
the system of non-linear difference equations, obtained by the 
classical five-point scheme at interior grid points and the 
elementary approximations (3.1.12) or (3.1.13) to the boundary 
equation, possesses solutions. Moreover, there exists a maximal 
and a minimal solution, u, v, respectively, such that any solu- 
tion z of (3.2.1) satisfies 

v < z < u. 

5. 5 Uniqueness of the Solution of the Difference Equations . 

Add the additional hypothesis that y is non- increasing 
in the dependent variable, i. e. , if 

w < z then y(z) < y(w) , (3.3.1) 

(2.2.14) being a special case. Suppose that u and _v are 
two different solutions of (3.2.1) defined by (3.2.24) and 
(3.2. 27) respectively. Define w by 

w = u - v. (3. 3. 2) 

Then at all points of 

w > 0 (3. 3. 3) 

where, for at least one point, w^ > 0 by assumption. In addi- 
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tion if this were not true for some point of P^g a contra- 
diction is immediate since in that case. 

Aw = 0. (3. 3. 4) 

Thus it may he assumed that strict inequality holds for at 
least one point on P^g. For such points Pq 

w 0 - w 2 

— a — = r( uq) - r(v 0 ) 

< 0 (3.3.5) 

hy ( 3. 3. 1) . So that 

Aw < 0 (3. 3. 6) 

This with (3.3.3) implies w = 0. That is, 

u = v 

and since any solution z_ of (3.2.1) must lie between them, . 
uniqueness is proven. 

Thus there is the result analogous to that of the con- 
tinuous case: 

Theorem 5. 5. 1 . - If the conditions for the existence of 
solutions of (3.2.1)' are supplemented hy 

r u <°, 

then there is only one solution of the non-linear difference 
equations. 

5. 4 The Practical Solution of the Difference Equations . 

Although section 3. 2 covers existence hy a proof which is 
constructive, if h is very small as compared with the diameter 
of R, the algorithm calls for the repeated direct solution of 
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large systems of linear equations. For one- dimensional prob- 
lems and, if h is not too small, for two-dimensional problems 
as well, the method is practical for use on large computers to 
obtain the solution of (3.2.1). 

Another possibility also requiring repeated direct solu- 
tions of linear systems is Newton's method set up for (3.2.1) 
analogously to the way it was done in section 23 for the con- 
tinuous problem, again under the assumption that y u < 0. The 
advantage of this method is its convergence properties, namely 

h - im+lH ^ K h - ^tnll 2 (3.4.1) 

whereas the first method's convergence is only linear. A dis- 
advantage of this latter method is the fact that the matrix of 
coefficients of the linear system of equations is different for 
each iteration. Therefore recourse cannot be made to "factor- 
ing methods" (see Varga [25]) to reduce the number of arith- 
metic operations per iteration as can be where the matrix is 
independent of m. 

A third possibility is suggested by the "simplified 
Newton's method" of Collatz [8]. Here (examining the continu- 
ous problem, the discrete case being entirely analogous) the 
factor T I (u m _2) appearing in (2.3.6) is replaced by y'(uo) 
and not changed from iteration to iteration as it is in the 
ordinary Newton's method. This third method appears "inter- 
mediate" to the other two in that the advantages of the first 
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in the invariance of the iteration matrix are retained and, 
while the convergence properties are not those of the ordinary 
Newton's method, they are surely better than those of the first 
method if Uq is already a reasonable approximation to u. 

Requiring solutions with a fair amount of detail to two- 
dimensional problems and with a moderate amount of detail to 
three-dimensional problems is becoming increasingly common in 
technological applications. In view of this the previously 
mentioned methods of obtaining such solutions for (3.2.1) no 
longer become practical. Instead it is necessary to make use of 
some sort of "relaxation technique" for these large order non- 
linear systems. 

There are some results available on the convergence of 
various relaxation methods for systems of non-linear algebraic 
(or transcendental, depending on y) equations. See, for ex- 
ample, Lieberstein [15], and Bers [5] (where here the systems 
arise from a difference method for Dirichlet's problem for 
Poisson's equation with non-linear source term) and 
Schechter [23]. 

It is shown here with the aid of Schechter' s results, how 
(3.2.1) may be solved by a (non-linear) "Gauss-Seidel" type re- 
laxation and that even "over- relaxation" is possible. In fact, 
the numerical results of section 7. 1 were obtained by using the 
one-line successive over- relaxation iterative method where only 
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one equation per mesh line, corresponding to the mesh point on 
the boundary Tg* was non-linear. Thus the tri-diagonal matrix, 
which couples unknowns along the line, is partially factored 
into a product of upper and lower triangular matrices. The 
factoring per iteration can then he completed by solving exactly 
this single non-linear equation in a sihgle variable. 

To apply the results of Schechter's paper directly, it be- 
comes necessary to take advantage of the smoothness of rg to 
transform R into another region R^ whose corresponding 
boundary is a subset of the x^-axis, Xg = 0. Upon making 

the transformation, the discretization is carried out by means 
of a square grid whose elements are parallel to the axes, so 
that the difference equations on Pg^ are all of the form 
(3.1.12). Having done this, it is possible, after rejecting 
the identity difference equations on P-j-p,, to make A a sym- 
metric matrix, which is the point of the above transformation 
and a requirement for the application of Schechter's work. 

It will be discussed in Chapter V, section 5.4 how sym- 
metry of A (truncated) may be obtained without insisting upon 
such a transformation. 

However symmetry is obtained, i. e. , whether by the above 
transformation, or by the use of other boundary approximations. 


let 
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G(u) = ^u^(A)u--i y ^ P r(P±,v)dv , (3.4.2) 

P l“ r 2h 0 

where here it has "been assumed that Tg is a mesh line and has 
been divided exactly into subintervals of length h. If the 
points of ^ H Tg turn out not to be nodes of the grid, the 
above summation is modified accordingly (being a trapezoidal 
approximation of a boundary integral) . 

A is positive definite. If B is defined by 


b.. = d 2 G _ 

1J dui duj 


“lj - T? r T (P k ,v) 


V- 


U k> 


since r u < 0, B is positive definite along with A. 

r(u) = grad G(u), 

Then (3.2.1) is equivalent to 


(3.4.3) 


Define 
(3.4. 4) 


r(u) = 0 , (3.4.5) 

It can be seen that for u e the real n-dimensional vector 
space, 

G(u) >4m,NX Z(r 2 ) (3.4.6) 

tr 


and fijj is a "solvent set". That is, [23], having chosen 
arbitrarily u e f^, u^ may be altered, determining a new 
u' s such that r^(u’) = 0. Here X(a) denotes eigenvalue 
of A and 2(Pg), the length of Pg. was defined prior 

to (2.1.25) and N by (2.1.30). 

The iterative procedure, analogous to the Gauss- Seidel 
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process for systems of linear equations, proven convergent in 
Schechter [23], is given as follows: 

Let u( m ) e iftn. Define u( m+ l) "by 



4 m+l) = k ^ (3.4.7) 

That is, the i^* 1 equation of (3.4.6) is regarded as a function 

of the i^* 1 variable and is solved for u^ with guessed values 

(previous iterate) for all other u.(j ^ i) in that equation. 

J 

Schechter actually proves that any "free steering" method of 
changing i converges, hut as a special case, the non-linear 
Gauss- Seidel method, consisting of taking each equation of 
(3.4.6) in turn, does converge. 

Schechter also proves a theorem on convergence of a type 
of over- relaxation, hut the satisfaction of the conditions of 
his 'ftheorem depends very strongly upon the x at hand and in- 
volves some greater restrictions. 

The results of this section may he summarized in, 

Theorem 3. 4. 1. - If the matrix A of coefficients of the 
system of finite difference equations (3.2.1) is symmetric and 
if r u is non- positive and satisfies the conditions of the 
existence theorem 2.1.2, then the solution of the non-linear 
system of equations may he determined hy a non-linear Gauss- 


Seidel iterative method. 


CHAPTER IV 


CONVERGENCE OF SOLUTIONS OF THE DISCRETE PROBLEMS 
TO THE CONTINUOUS 

That the solutions of the Boundary value problems (2.1.1 - 
2. 1. 3) may he approximated arbitrarily closely by the solutions 
of the non-linear system of equations (3.2.1) is proven first 
under the restriction y < 0. It is then shown that the ele- 
ments Ujjj of the sequence generated by (2.1.33 - 2.1.35), 
(2.1.40 - 2.1.42) can be approximated as closely as desired by 
the corresponding elements of the sequence generated by 

(3.2.3), (3.2.7). Finally convergence of u of (3.2.24) to u 
of (2.1.51) as h -► 0 is shown under assumption of uniqueness 
of the solution. 

Actually, to show convergence as h goes to zero, the 
approximating system is modified by rejecting equations written 
at points in a circular neighborhood of the corners where the 
"Dirichlet boundary" and the "Neumann boundary" intersect. The 
radius, k, is taken larger than h but goes to zero with h. 
This approach is similar to the one used by Roudebush [22] and 
Rosenbloom [21] where equations were rejected near all of a 
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non-smooth boundary (assumed merely to have strong barrier 
function) to show (see [22], p. 115) that the discrete approxi- 
mation approaches the solution, assumed Holder continuous with 
exponent a, of a Dirichlet problem for a uniformly elliptic 
differential equation, at the rate h a2 / 2 ( 1+a ) which is h 1 / 4 
if u satisfies a Lipshitz condition. 

This modification imposes no hardship on the results of 
the other chapters. 

4. 1 Convergence When y u < 0. 

It is now assumed that t u < 0, that all the smoothness 
hypotheses (2.14 - 2.1.9) on y are satisfied and that 
and Tg are both smooth, the condition on r-j_ being later 
relaxed to get slightly weaker results. 

We choose h and k so that 


0 < h < k, (4. 1.1) 

both tending to zero in what follows. A modified grid is now 
defined where some of the grid points near the corners O r 2 
are dropped from consideration, by setting 


R$ = P e Rv 


d(P,r 1 n r 2 ) > kj 




P k = 


fP e R h - Rg 


there is a 

Q € Rg 7 

such that | 

Q - P| < hj 


R h, k " R k U r k" 


(4.1.2) 


For these definitions it is helpful to examine figure 4 which 
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is an illustration relative to figure 1 of R^ k where k has 
"been taken equal to 2h as an example. The sets R h?k , r h,k> 

-)f 

r h,k^ r h,k, 1> etc * have -fche same relation to R h ^ k that R h , 



T hfl , etc. respectively, have to R k , adding the points 
of r£ to the set r k k 

A new finite difference houndary problem is defined, in 
R^ k by setting 

UkCP) =0 if P e T k (4.1.3) 

retaining the original equations at all other points of R k k< 
Its solution is guaranteed unique by Theorem 3. 3. 1. 

The difference function Z is now formed at points P^ of 

\,k : 

Z = - u (4. 1. 4) 

where u^ = u(Pj_) and u is the unique solution of (2.1.1 - 
2. 1. 3) guaranteed by the results of CHAPTER II. To estimate the 
error in the part of the matrix operator A of (3.2.1) 

which still remains will be applied to Z. For later purposes 
call the "truncated" matrix A^ Thus from (3.2.1) 

^k£ = ^kUk " A kii 

= dj,(U k ) - Ak u. (4.1.5) 

Let D 2 v(Pj_) denote any second partial derivative of v 

taken at the point P^. Then results of Miranda [17] show that 

there is a constant M 2 such that 

|d 2 u(P;l) - D 2 u(P 2 )| < — § | P x - P 2 | a . 

k^ 


(4.1.6) 
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To examine the error Z at a point P Q e R h)k U T* k by 
considering the "formal error" A^Z, figure 2 should be con- 
sulted. A Taylor's series procedure similar to that of 


Roudebush [22] is used: 

u l = u O + h l u xO + i h l^xxl C 4 * 1 - 7 ) 

u 2 = U 0 + h 2 Uy 0 + \ h 2 ^2 (4.1.8) 

u 3 = u 0 - 1^0 + | h 3 Scx3 (4.1.9) 

u 4 = u 0 - VyO + | V (4.1.10) 


where, for example, u^ is c)u/5x evaluated at the point Pq 
and u xxl is S 2 u/dx 2 at seme point between Pq and P^. 
Subtracting Uq from both sides of (4.1.7) and (4.1.9), divid- 
ing (4.1.7) by h]_ and (4.1.9) by h^, and adding one obtains 

U 1 ~ u 0 , u 5 - u 0 _ h l ~ . h 3 ~ _ ... 

hj_ hj 2 u xxl 2 u xx3 * (4.1.11) 

Multiplying by 2/(h^ + hj) and subtracting u xx q one obtains 

2 / U 1 - U Q u 5 - up \ h i . 

^ + h 3 y hj_ + h-z; J ~ ^xO “ h + ^ ' U xxl “ u xxo' 

+ +"h 3 ^^xx3 - u xxO^ (4.1.12) 

Analogously, operating on (4.1.8) and (4.1.10) one obtains 
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2 f U 2 " U 0 + u 4 - U Q^ 


h 2 + h 4 


h«- 


Ci 

'Vyo = h + h ‘ ^ u yy2 " 'Vyo) 


2 4 


h 4 

+ h 2 + h 4 ^^4 ” ^yyo) (4.1.13) 


Adding (4.1.12) and (4.1.13) one obtains "by (4.1.6) 

M 


(A k u)(P 0 ) +V 2 u(P 0 )|< 


_2 

.2 


hp _ h 


, h l + h 3 


a . "3 


p i- p or + i^-Thj 


P 3 “ P o| 


a 


Sfcfc- 


2 4 


■0 


a + h 4 


h 2 + h 4 


P 4 - P 0 


a 


s%*. 

k 2 


(4.1.14) 


To examine the formal error at points of figure 3 

should "be consulted. 

Let u* he the value of u where the normal intersects 
the grid line between the points labelled 1 and 2. Then 


u~ = 


= u* + d 

SnJ , 


d£ ^u 

2 3n 2 <1o) 


as before. Also 


u 2 = u* - d 2 


up = u* + d x 



u 

2 5y 2 



d i a 2 u , . 

+ -s? Ui) - 


( 4. 1. 15) 


(4.1.16) 


(4.1.17) 


Eliminating (du*/dy) between these equations yields 
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diu 2 + d 2 u x 

. d l d 2 

a 


4- rl 


°* d l + d 2 

- 2(d 1 + d 2 ) 

d 2 

ay 2 (na) 

+ d x 


. d l d 2 M2 

s_s_ 7 


(4.1.18) 


so that 


By (4.1.15), 


’ll-#- cl^u.*^ 

d d(d^ + dg) 


d l d 2 M2 
- 2d k 2 


d 2 m 2 


M = !° - Sul + 0 /a\ 

d d H 0 V/ 


(4.1.19) 


(4.1.20) 


where d < V2 h. Combining (4.1.19) and (4.1.20) one obtains 


for P 0 e r h,k, 2 ’ 

(A k u)(P 0 ) = r(u 0 ) + o(±y 

Therefore 

(%z)(p 0 ) = r[u k (P 0 )] - r(u 0 ) + o(±) 

= r’(u)[u k (P 0 ) - uq] + o(A), 

i. e. 

(A k z)(p 0 )| = o(A), 


(4.1. 21) 


(4. 1. 23) 


(4. 1. 24) 


where here has a slightly different meaning than it has in 

section 3. 2 in that the diagonal elements affected have positive 
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numbers added which are not constant. Nevertheless, the in- 
equality (3.2.11) still holds, and this is the key to the error 
estimate in the non-linear case where y u < 0. A similar re- 
sult holds for the simpler boundary equation (3.1.12) wherever 
it may be applied. Here (4.1.17), with the subscript re- 
placed by "0" leads to a result like (4.1.24). 

For Pj_ e rg it is already true that 

|Ufc(Pi) - u(P i )| = 0( k a ) . (4.1.25) 

For Pi e T h)k>1 - rg 

U k (Pi) = u(P ± ) = 0. (4.1.26) 

Compiling (4.1.14), (4.1.24), and (4.1.25) yields 

Aj£ = 0(k a ) + 0^^ (4.1.27) 

as h-»0. Letting k = hP and maximizing the order of con- 
vergence in (4.1.27) yields finally |3 = a/(2 + a) and thus 
there is a constant M > 0 such that 

| A^Zj < Mh a2 / ( 2+a ) e (4.1. 28) 

where e is a vector, all of whose components are 1, 

and |x| is defined by 
Since both A^ 1 and A^ 1 


e = 


(4.1. 29) 


X = I Xj 


(4.1.30) 


are positive matrices, (4.1.28) and 
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(3.2.11) imply 

J Z| < Mh a 2 /(2+a)^-l £ 

< Mh a2 / ( 2+a ) Aj^e , ( 4. 1. 31) 

To complete the proof of convergence, it must he shown that the 

-1 

components of e are hounded independently of h. 

Let £ he the unique solution of (2.1.8 - 2.1.10), guar- 
anteed to have the same smoothness as u as a special case of 
Miranda's theorem. 

The vector £ is defined hy 

i ± = £(Pi). (4.1.32) 

For Pq e R k k , applying (2.1.12), one obtains, in the manner 
that (4.1.14) was found, that 

(\I)( p o) = 1 + 0 [h a2 / ( 2+a ) , (4.1.33) 

so that, for all sufficiently small h, one has 

(%£H P 0 ) >!• (4.1.34) 

Also, for Pq e k 2 , applying (2.1.14), one obtains, in the 
manner that ( 4. 1. 24) was found, that 

(A^)(P 0 ) = 1 + 0 jh a / ( 2+a ) (4.1.35) 

so that for all sufficiently small h, one again has 

(4d)( p 0 ) (4.1.36) 

This is surely true for Pq e Ph,k,l since here, A^ is the 
identity. Therefore, 
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Aj^-e < 2£, 

and by ( 4. 1. 31) , 


(4.1.37) 


|Z| < 2Mh a2 /( 2+a )^. (4.1.38) 

Now, since £ is a bounded function, we have thus proved: 

Theorem 4.1.1 . - Under the hypotheses of the uniqueness 
theorem. Theorem 2.2.1, if k = h a / ( 2+cc) f converges to the 

solution u of (2.1. 1 - 2.1.3) at the rate 



In particular, if r u and the second derivative with respect to 
are parameter of the representations of and r 2 satisfy 

Lipschitz conditions (a - l) with respect to arc parameter, 
then 

|U k - uj = Oft 1 / 3 ). 

4. 2 Convergence of the Discrete Iterate to the Corresponding 
Continuous Iterate . 

Let Uq and Vq he respectively defined by (2.1.33 - 
2.1.35) and (2.1.55 - 2.1.57) with N' = N + 1. The sequences 
{u^} and (v n ) are then defined by (2.1.40 - 2.1.42), starting 
with Uq and Vq respectively. Let Uq and Vq be the 

solutions of (3.2.3) and (3.2.25) with % replaced by Rjyp 
which implies the replacement of A, b and b by Ap, bp, 
and bp, respectively. The sequences {U^} and { V,^} are then 
defined by (3.2.7), again with A, bQ^) , and b( V m ) replaced 
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by , and b v (V m ) , respectively. 

The results of section (4.1) already show that u Q and 
Vq may he approximated arbitrarily closely by Uq and Vq. 

It remains to show that this is also true for ^ and v n for 
any n. 

Assume that u n _^ and v n _j_ may be approximated by 
U ri _^, and V ri _i , respectively, as closely as desired. It will 
be proven that u^ can then be approximated arbitrarily closely 
by U n , the argument for v^ being completely analogous. Thus, 
by hypothesis, for arbitrary e > 0 there is an h^ < 0 such that 
for all h satisfying 0 < h < hQ 

IHn-1 - “n-ll < ( 4 - 2 -D 
where the subscript, k, appearing in the previous section, has 
been dropped for simplicity and will not appear modifying func- 
tions or sets. Let W n be defined by 

A In = £(u n _i,W n ) (4.2.3) 

J, T_ 

where the i u component of b is zero unless Pj_ e Here 

i = r(P i ,u n _ 1 ) + k|w n (P.) - u^U^. (4. 2. 4) 

But by the results of section 4. 1, there is an h^ > 0 such 
that for all h satisfying 0 < h < h^_ 

|W n - u n | < €. (4.2.5) 

Let 




(4. 2. 6) 
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Then 

A k£ = A kU - AkW n 

= b(U n ) - 

= i(V^ n ). (4.2.7) 

where b^ = 0 unless P^ e r^g. Here 

^V^Hn) = rCPi^.i) + k[Un( p i) - U n -l( p i^ “ 

- k ^ n (Pi) - Un.iCPj)]. (4.2.8)' 

Thus at these points, 

|(A k fi)(P i ) - k%| < ItKp^u^) - 

+ k |Un-l( P i) - "n-l( p i)| 

< 2k|U n _ 1 (P i ) - u n _ 1 (P i )|. (4.2.9) 

Therefore if h is such that 0 < h < min(hQ, h^), then 

| Aj£| < 2k e e (4.2.10) 

so that, as in section 4.1, 

|%| < 4k € ^(Pi). (4.2.11) 

Recalling that £ is a hounded function, the following theorem 
holds: 

Theorem 4. 2. 1 . - Without the assumption T u < 0, the iter- 
ates u n for the continuous problem (2.1.1 - 2.1.3) and the 
iterates n for the discrete problem A^U^, = ^(U^) are 
considered. If k = h a / (2+u) , then 

ia k , n - %l - op/ ( H, 

for n > 0. Similarly, this is true for the iterates V^, n and 
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4.5 Convergence of the Discrete Limit Function to the Continuous 
Limit Function 

Under the assumption that the solution to (2.1.1 - 2.1.3) 
is unique (hut not necessarily assuming that r u < 0) it is now 
shown that, if is defined hy = d k ( U- k .) where k = 

h a /(2 + a) and u ^y (2.1.51), u can he approximated as closely 
as desired hy U- K hy letting h go to zero. 

If v is defined hy (2.1.58), then uniqueness of u 
implies that 

u = v (4.3.1) 

Although not done in the previous sections (e.g. , 4.2, 3.2) 
it now becomes necessary to subscript U and U n with k. 

That is, in CHAPTER III, where they were first discussed, 
and U were defined for a particular h. Let e > 0. Then 
from the results of CHAPTER II prior to Theorem 2. 1. 2, there is 
an Nj such that for all P e R^^. and n > N, 

u < u n < u + -|. (4.3.2) 

Ey section 4. 2 there is an h-j_ > 0 such that for all h with 
0 < h < h^ and all P.^ e 

l«n,k-“nl < h < 4 - 3 - 3 > 

Therefore, for these n and h 

« k < u„,k < % + | S u + «• 


(4.3.4) 
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Similarly there are N 2 , hg such that n > N 2 , 0 < h < h 2 
implies 

( 4 . 3 . 5 ) 

Equation (4.3.4) implies 

Ufc - u < e. (4.3.6) 

Equations (4.3.5) and (4.3.1) imply 

u - Ik = v - Sc < v - 3k < e * ( 4 - 3 - 7 ) 

But (4.3.6) and (4.3.7) are the results desired. That is, if h 
is such that 0 < h < min(h^,hg), then 

|Uk ‘ u l < € (4.3.8) 

at all points of k - 
The result is 

Theorem 4. 5. 1 . - If there is only one solution, u, of 
(2.1.1 - 2.1.3) then converges to u as h -*■ 0 if 

k = h a /( 2+a ), i. e. , for points of R^ k , 

|U k - u| = o[h a2 /( 2+a )] . 



CHAPTER V 


GENERALIZATIONS 

The results in this thesis thus far are stated for the 
Laplace equation in two dimensions for mixed boundary condi- 
tions which are in part Dirichlet. Extensions to more general 
elliptic equations, with non-linear source terms, extensions 
to higher dimensions, and extensions to more general assumptions 
on the boundaries and boundary conditions can be carried out in 
some spots very simply in others with some difficulty. 

Although the numerical analysis was carried out on the basis 
of the classical five-point difference approximation to the 
partial differential equation on square networks, it can be 
extended in the various directions taken in recent years by 
workers (for example Bramble and Hubbard [7] and Roudebush [22]) 
in this field. The purpose of this chapter is to indicate the 
directions of generalization. 

5.1 Higher Dimensions . 

Although the discussion in this thesis was carried out 
assuming everywhere that the equations were in two independent 
variables, almost all arguments could have been carried out for 
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differential equations in n-dimensional space. One reason for 
not writing it all for n variables is the greater length of 
the expressions which would result in the numerical treatment 
starting in CHAPTER III, as, for example, for (3.1.9) and the 
correspondingly larger increase in detail in the arguments for 
the error estimates, as, for example, the arguments leading 
from (4.1.7) to (4.1.14). Most of the references are either 
already in n-dimensions, have statements indicating the exten- 
sion of their results to n-dimensions, or can be extended to 
n-dimensions . 

As an example, consider Miranda's work [17] which is basic 
to this paper. The results here are stated in m-dimensional 
space and carried out in great detail under the assumption 
m > 3, with statements about where arguments have to be modified 
for the case m = 2. Two reasons for this are the difference 
in form of Green's functions appearing in them, and the neces- 
sity of considering the hypergeometric function for m > 3, 
the case m = 2 being much simpler. 

Probably the greatest difficulty would arise in the deriva- 
tion and treatment in the error estimates of the difference 
equations on p„ where the non-linear boundary conditions 

C 

apply, if they are to be set up analogously with the two- 
dimensional treatment here. 

5.2 Partial Differential Equations . 


Although Laplace ' s equation was treated throughout this 
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paper, immediate extensions can be made to more general homo- 
geneous equations and to their non-homogeneous counterparts. 

In fact the paper upon which this thesis rests most heavily, 
Miranda [17], treats partial differential equations of the form 
m m 

^ , , a 2 

a ik (x) 

i,k=.l i=l 

where (a^) is positive definite in all of R and satisfies 
conditions similar to those of uniform ellipticity. The coef- 
ficients in addition satisfy certain Holder conditions in do- 
mains not containing O Tg permitting certain orders of 
approach to infinity as the distance of these domains from 
r O r 2 approach zero. 

Again, the greatest difficulty appears to be in the numeri- 
cal analysis and will require some relaxing of the class of 
difference approximations . This aspect is treated in section 
5.4. 

5.5 Boundary Conditions . 

Miranda indicates some ways to weaken the smoothness re- 
quired of cp (and thus for the $ m of section 2.l) in the 
boundary conditions (corresponding to (5.2.l)) 

w*) tr + = <pW 

K. 

k=l 




+ °( x ) u = (s- 2 - 1 ) 


(5.3.1) 
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on Tg • Pursuing this would lead, to convergence (in the sense 
of CHAPTER IV) of a lower order. 

Stampacchia [24] discusses various weak ways of satisfying 
the boundary conditions and discusses non-linearities in both 
partial differential equation and boundary conditions from the 
point of view of calculus of variations . 

With regard to Fq there are also some directions of 
generalization : 

(i) The type of boundary condition may be other than 
Dirichlet. For example 

= cc(u Q - u), a > 0, (5.3.2) 

is a possibility and does not lead to too much difficulty in 
the numerical analysis as monotone properties are still present 
and boundedness of the inverse independently of h still holds. 
This is an especially useful generalization in applications to 
heat transfer. 

(ii) P]_ = 0 . In other words, the non-linear Neumann 
type of conditions is permitted all around. Here some addi- 
tional restrictions need to be made on y for existence as 
well as uniqueness. 

(iii) Relaxation of smoothness of Tq : Here the results 

of this paper need to be combined with those of Roudebush [22] 
and results obtained are the essentially 0(h^-/^) convergence 


obtained there . 
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5.4 Difference Equations . 

Here, a generalization can quickly be made since it will 
be recalled that (3.1.11) was not even analyzed in CHAPTER IV. 
That is, the mesh might as well have been assumed not neces- 
sarily square, but rectangular, to begin with, gaining the 
advantage of possibly "fitting" the domain R a little better 
by R^. Convergence is then analyzed on the basis of 

h = max hj_. (5.4.1) 

Other mesh patterns have been considered in the past 
with varying degrees of success . It might be advantageous 
to consider them here since an arbitrary triangulation of R 
could possibly "fit" it better and is amenable to differencing 
based upon the integration technique. (See Varga [26], p. 184). 

It would certainly be advantageous to consider the inte- 
gration technique since it takes boundary conditions with 
normal derivatives without difficulty and produces automatically 
a symmetric matrix. This is the ideal situation for the appli- 
cation of Schechter's iterative method's. Some error analysis 
needs to be carried out yet for this type of boundary condi- 


tion treatment. 



CHAPTER VI 


SUMMARY AND CONCLUSIONS 

This thesis is concerned with solving Laplace's equation 
with mixed boundary conditions, part Dirichlet and part Neumann, 
in bounded two-dimensional regions. On that part of the 
boundary where the normal derivative enters, it is a non-linear 
function y of the solution itself. Such boundary value 
problems have solutions under some, not too restrictive hypoth- 
eses on the smoothness of y, including the condition 

lim sup Ij&ul <Q (e.i.l) 

and on the smoothness of the boundary. 

A constructive existence proof is carried out giving two 
solutions, (possibly equal) the largest possible and the 
smallest possible, as limits of two sequences of functions 
satisfying two sequences of linear boundary value problems 
related to the original. When 

T u < 0, (6.1.2) 

there is only one solution and thus the two limits are equal. 

The intention of this research was to prove that the 
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application of finite differences could be made to elliptic 
equations with non-linear boundary conditions, with a background 
of mathematical certainty that the functions so obtained are 
near those which are the solutions of these boundary value 
problems. This is carried out rigorously in the fourth chapter, 
using the most elementary differencing method based upon the 
classical five point scheme, as modified in the third chapter 
to fit problems of this type. Although the exponent is not 
large, suggesting the need for large numbers of points to get 
moderate detail in the solutions, numerical results suggest 
that with more care, better convergence of the solutions of the 
difference equations to that of the original problem is obtained 


in practice. 



CHAPTER VII 


NUMERICAL EXPERIMENTS 


7 .1 Radiation of Heat to Space . 

The following example arises from the necessity of rejectr 
ing waste heat; into space. It represents a section of a rec- 
tangular radiating fin. 

It is required to find approximations to the temperature 
field, considered essentially to be two-dimensional and satis- 
fying the following equations 


V 2 T = 0 

in 

R: 0 < x, y < 1, 

(7.1.1) 

ST A 
= 0 on 
on 

r 0 : 

x = 0, 0 < y < 1, 
y = 0, 0 < x < 1, 

(7.1.2) 

T = 1 on 

i — 1 

Ch 

x = 1, 0 < y < 1, 

(7.1.3) 

f-T 4 , T > 

' 0, T < 0 

on r 2 : y = 1, 0 < x < 1, 

(7.1.4) 


ST 

ST 


in what is really only the first quadrant portion of a square 
two units on a side, but considered as above because of the 
symmetry of the original problem. 

Define 


C 1 (x,y) = 2 - x 2 + | y 2 . 


(7.1.5) 
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which satisfies 


V 2 ^ = -1 in R, 

(7.1.6) 

asi 

dn ~ 0 on r 0' 

(7.1.7) 

= 1 + | y 2 on r x , 

(7.1.8) 

*S, 

sir " 1 on r 2 • 

(7.1.9) 


Therefore, if £ is to he defined by (2.1.12 - 2.1.14) on the 
entire 4-quadrant configuration, (7.1.6 - 7.1.9) implies that 

5(x>y) < ? x (x,y) (7.1.10) 

for (x,y)eR. 

To place this in a form where the a priori estimate of 
section 2.1 may be obtained, it is necessary to define T^ by 

T = T x + 1, (7.1.11) 

and 

r -(T x + l) 4 T 1 > -1 

r(T 1 ) = < (7.1.12) 

0 T x < -1. 

Problem (7.1.1 - 7.1.4) is then rewritten as 

V 2 ^ = 0 in R, (7.1.13) 

ST 

ZT = 0 on r 0 , (7.1.14) 


T x = 0 on r x , 


(7.1.15) 
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dT 


= r(T x ) on r 2 . 


(7.1.16) 


Then for -1 < T x < 0, 

T(T, ) 


? ■ ■ ' 'V • *'• (- *) 


- l T ll ■ 


Therefore, if we take 


n 

|u| 


l u j < 1 


Tj(u) = / 


(7.1.17) 


'TTp M >X > 


where m > 1 is free to be chosen, all conditions of section 
2 .1 are met. 

Since 

max £ = 2.5, (7.1.18) 

r 2 

it is possible to take M-j_ any number such that 

Mj_ > 1, (7.1.19) 

simply by making m large enough in (7.1.17). Thus 

M 3 = 2.5, (7.1.20) 

i.e. , 

| T !j < 2.5, (7.1.21) 

so that 

|T| < 3.5. (7.1.22) 
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This is a coarse estimate in this case, which is subject 
to some improvement, perhaps, if £ were used instead of 
However, applying the strong maximum principles (see Bers, 

John, Schechter [5], p. 151) shows that 

| T| <1 in R. (7.1.23) 

The boundary value problem was replaced by a system of 
finite difference equations following the methods of section 
3.1. A square mesh of mesh width h was used and was so 
situated in the square that the edges were parallel to the 

•ft 

grid lines and of distance h from the resulting 1^. 

Equations near the corners were not rejected as required 
for the convergence proof of CHAPTER IV but left in for ease 
of programming. The numerical results presented show that, at 
least in this case, convergence as h -» 0 is still possible. 


h 

Error 

Apparent 

order 


0.01211 



l/ 8 

.00797 

1.52 

l/l6 

.00468 

1.70 

l/32 

.00264 

1.77 


The third column of this table is obtained by replacing (by 

assumption) the inequality implied in (4.1.39) by equality 

and the exponent by an unknown, a the "apparent order", 

8 

namely we assume 

n n _ a a 

||error|| = e = Mh 

and solve for a , based on the second column. 

8 
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As already announced in section 3.4 this example was 
solved by a line-by-line successive over-relaxation scheme 
(S.O.R.). A series of experiments on some rectangular problems 
indicate that, as in the linear case where optimum values are 
obtainable rigorously, there is a value for o>, the S.O.R. 
parameter, which is optimum. This is taken in the sense that 
any other value leads to more iterations before the error in 
an aribtrary initial approximation is reduced to an acceptable 
level. 

In conclusion, this numerical example has strengthened 
this thesis in indicating, that the predicted convergence of 
approximate solutions to the solution of the continuous, does 
indeed take place in practice. However, it serves to point 
up its weaknesses and thus indicate further directions for 
research. First, it should be possible to prove that one need 
not reject equations near the corners and an approach using 
"smoothing" operators might apply here as it did in Parter E 18 ] . 
Secondly, something should be possible in analytically pre- 
dicting a value which is optimum for the relaxation parameter. 
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